On the Scission Point Configuration of Fisioning Nuclei 
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The scission of a nucleus into two fragments is at present the least understood part of the fission 
process, though the most important for the formation of the observables. To investigate the potential 
energy landscape at the largest possible deformations, i.e. at the scission point (line, hypersurface), 
the Strutinsky's optimal shape approach is applied. 

For the accurate description of the mass-asymmetric nuclear shape at the scission point, it turned 
out necessary to construct an interpolation between the two sets of constraints for the elongation 
and mass asymmetry which are applied successfully at small deformations (quadrupole and octupole 
moments) and for separated fragments (the distance between the centers of mass and the difference 
of fragments masses). In addition, a constraint on the neck radius was added, what makes it possible 
to introduce the so called super-short and super-long shapes at the scission point and to consider 
the contributions to the observable data from different fission modes. The calculated results for 
the mass distribution of the fission fragment and the Coulomb repulsion energy "immediately after 
scission" are in a reasonable agreement with experimental data. 

PACS numbers: 02.60.Lj, 02.70.Bf, 21.60.-n, 21.60.Ev, 25.85.Ec 

Keywords: nuclear fission, deformation energy, scission point, super-long, super-short shapes, mass distribu- 
tion, total kinetic energy 
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I. INTRODUCTION 

The shape of a nuclear surface is a basic notion in many 
theoretical models of nuclear structure and reactions. A 
good choice of the shape degrees of freedom reduces sub- 
stantially the computation time and is often, especially 
for the description of fission process or fusion-fission re- 
actions, a key to the success of the theory. 

In past a lot of shape parameterizations were proposed 
and used. One class of shapes relies on the expansion in 
a complete set of functions like the expansion of the ra- 
dius vector R(9) [l| or the profile function squared p 2 (z 



in Legendre polynomials. In the parametrization 
the deviation of the shape from the basic Cassini ovals is 
also expanded in Legendre polynomials. Another possi- 
bility is given by the introduction of a restricted number 
of deformation parameters, like in the parametrization 
of three smoothly joined quadratic surfaces Q, the two 
center shell model d, the Funny-Hills parametrization 
Q or modified Funny-Hills parametrization Q. 

All these shape parametrizations are restricted to a cer- 
tain class of shapes. In all these cases the question arises 
whether the given class of shapes is complete enough to 
represent the essential properties of the investigated pro- 
cess. 

A method to introduce the shape of the nuclear sur- 
face which does not rely on any shape parametrization 
was suggested by V.Strutinsky already in In this 

method one defines the profile function p(z) of an axi- 
ally symmetric nucleus by the minimization of the liquid 
drop energy with respect to the variation of p{z) under 
additional constraints which fix the volume and elonga- 
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tion of the drop. However, due to numerical difficulties 
this method was not widely used in the past. 

Only recently [l(| it turns out possibly to solve the 
variational problem of Q in a very broad region of de- 
formations ranging from a disk (even with a central de- 
pression) to two touching spheres. The fission barriers 
calculated by this method [ll| were found to be in a rea- 
sonable agreement with the experimental results. 

In the present work the potential energy landscape is 
investigated for larger deformations - at the scission point 
(line, hypersurface). The scission of a nucleus into two 
fragments is at present the least understood part of the 
fission process, though the most important for the for- 
mation of the observable data. 

For the accurate description of the mass-asymmetric 
nuclear shape at the scission point it turned out necessary 
to construct an interpolation between the two sets of con- 
straints for the elongation and mass asymmetry which are 
applied successfully at small deformations (quadrupole 
and octupole moments) and for separated fragments (the 
distance between centers of mass and the difference of 
fragments masses). In addition, a constraint on the neck 
radius was added, what makes it possible to introduce 
the so called super-short and super-long shapes at the 
scission point and consider the contributions to the ob- 
servable data from different fission modes. 

The paper is organized as follows. Section [TT] is a short 
overview of the Strutinsky optimal shapes prescription. 
The mass-asymmetric shapes are introduced in Sect. Mil 
The scission shape and the shape of separated fragments 
"immediately after scission" are defined in Sect. HVlTVl 
In Sect. I VII the super-short and super- long shapes are in- 
troduced and the calculated Coulomb repulsion energy of 
the fragments "immediately after scission" is compared 
with the experimental total kinetic energy of fission frag- 
ments of 236 U. Sect. IVIII contains a short summary. 
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II. THE OPTIMAL SHAPES OF FISSIONING 
NUCLEI 

The shape of an axially symmetric nucleus can be de- 
fined by rotation of some profile function p(z) around 
the z-axis. It was suggested in Q to define the profile 
function looking for the minimum of the liquid-drop en- 
ergy, Eld = E sur f + Ec ou \ , under the constraint that the 
volume V and the elongation R12 are fixed, 



MV -\2R12) = 0, 



with 



2tt 

V 



p (z)\z\dz , V = it I p (z)dz 



(1) 



(2) 



In (JTJ) Ai and A2 are the corresponding Lagrange multi- 
pliers. The elongation parameter R12 was chosen by [1] 
to be the distance between the centers of mass of the left 
and right parts of the nucleus, 

The minimization of .Eld — AiV ~ A2-R12 with respect 
to the profile function p{z) leads to an integro-differential 
equation for p(z) 



pp" = 1 + 00 -p[A 1 +A 2 |z|-f0x LD $s][l + (p' 



/\2i 



(3) 



Here i>s = &(z,p(z)) is the Coulomb potential at the 
nuclear surface, and ccld is the fissility parameter of the 
liquid drop [l2j], 



E (o) 

— -^Coul 



IE. 



(o) 

surf 



3 Z 2 e 2 
10 AttRIo 



49,4 



7T7T> ( 4 ) 



where a is the surface tension coefficient. In © and 
everywhere below the index refers to the spherical 
shape. 
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FIG. 1: The solutions of Eq. corresponding to different 
values of Lagrange multiplier A2 which fixes the elongation 
R12 



By solving Eq. © one obtains the profile function p(z) 
for given x^d and A2 (Ai is fixed by the volume conser- 
vation condition). The liquid drop deformation energy 

E^Jf = Eld — E^ (in units of the surface energy for a 
spherical shape) 



E A 



, def = E^/E^t = B sull - 1 + 2x LD (B Coul - 1) , (5) 
calculated for the shapes shown in Fig.[T] is presented in 

Fig.H In © Seoul = E C oul/4oul' B ^ri = E SUI {/E^ { . 
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FIG 2: Liquid-drop deformation energy ([5j as a function of 
the parameter R12 for different fissility parameters ild and 
the corresponding nuclear shapes at scission. 

One can see from Fig.[5]that the elongation R\2 of the 
shapes shown in these figures is limited by some maximal 
value R s \2- Above this deformation mono- nuclear shapes 
do not exist. This critical deformation was interpreted in 
Q as the scission point. Note that, at scission the neck 
radius is still rather large: the neck radius at the critical 
deformation is approximately equal to (0.25 — 0.30)i?o for 
a fissility parameter in the range 0.4 < xld < 0.9 

Another peculiarity of Fig. [5] is the upper branch of 
the deformation energy at large deformation. Along this 
branch the neck of the drop becomes smaller and smaller 
until the shape turns into two touching spheres. Both 
branches are solutions of Eq.©. It turns out, that the 
upper branch of E^ corresponds not to the minimum 
but to the maximum of the energy. Thus, it represents 
the ridge of the potential energy surface between the fis- 
sion and fusion valleys. 



III. THE MASS-ASYMMETRIC SHAPES 

The optimal shape approach of [8| can be generalized 
to mass-asymmetric shapes. For this aim one has to in- 
clude into Eq.© one more constraint fixing the mass 
asymmetry S of the drop, 

£ 

-(E LD -XiV- X2R12 - A 3 - . (6) 



Sp' 
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The mass asymmetry S is commonly denned by the dif- 
ference of masses Ml and Mr to the left and right of 
some point z*, 



Mr, — M 



R 



Ml + M r 



y I Sign(z 



z*)p 2 (z)dz. (7) 



In case that the drop has a neck, z* coincides with the 
position of the neck, z* = z n . By z n we mean here the 
point where p(z) has a minimum. For the pear-like shape 
the neck does not exist and z* could be defined in a dif- 
ferent way, see 10] . In the present work we are interested 
in the scission point configuration for which the neck is 
well defined. Then the Euler-Lagrange equation for the 
variational problem (|6]) has the form 

PP " = l + {p' f - p{XiV + X 2 \z-z n \ (8) 
+ A 3 Sign(z - z n ) + 10x LD * s (z)}[l + (p') 2 ? /2 ■ 

Eq.JHJ) can be solved in the same way as Eq.©. Some 

examples of the shapes at the scission point i?^ (max- 
imal possible value of R12) for few values of the mass 
asymmetry are shown in Fig. [3] 
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FIG. 3: The solutions of Eq. ([8} at the maximal elongation 
for a few values of the mass asymmetry 5 = 0, 0.1, ...0.5 
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The advantage of the variational problem in the form 
is that the constraints for the elongation and the mass 
asymmetry have a clear physical meaning. These are the 
distance R12 between the centers of mass of right and left 
part of the drop and the mass asymmetry 5 of the drop. 

The disadvantage is that due to the simplicity of the 
restrictions on R12 and 6, equation J5J contains some un- 
physical effects. Namely, because Sign(z — z n ) is a discon- 
tinuous function of z, the second order derivative p"{z), 
defined by Eq.©, and, consequently, the curvature of 
the surface, is discontinuous at z — z n , what should not 
take place for a liquid drop. Besides, for a pear-like shape 
the neck does not exists and it is not so clear how one 
could define z* in this case. Some possibility to define z* 



as the place of the largest curvature of the surface was 
suggested in [Toj . 

Besides R12 and S one could try another popular pair 
of constraints which are often used in the constrained 
Hartee-Fock calculations, namely, the quadrupole and oc- 
tupole moments, 

Q2 = yf dVr 2 P 2 (co S 6),Q 3 = i J dVr 3 P 3 (cosd) . 

(9) 

The moments Q2 and Q3 can be defined independently 
of whether the neck exists or not. The use of quadrupole 
and octupole moments as constraints 



— (£ LD - XiV - X2Q2 - A 3 Q 3 ) = 
dp 

leads to the following Euler-Lagrange equation 
pp" = l + {p'f-p%{z)[l + {p'f]^, 



with 



U{z) = X 1 V + X 2 



(10) 



(11) 



(12) 



\%{z — z cm ) 



( 2* Zcm ) 



2 3p 2 (z) 



WXLD^S- 



Eq. fTTl) can be solved in the same way as Eq.©. The 
examples of the scission point shapes for few value of the 
mass asymmetry are shown in Fig. |4j 
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FIG. 4: The solutions of Eq. (|11[) at the maximal elongation 
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for few values of the mass asymmetry S — 0, 0.1, ...0.8 



Comparing the shapes at the scission point calculated 
by Eq. (J5J and Eqs. dTT1)- p2l for the same mass- 
asymmetry one can see that these two sets of shapes are 
rather different. The Q21Q3 restrictions lead to scission 
shapes which are considerably "shorter" as compared 
with the scission shapes defined with i?i2,<5 restrictions. 
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The comparison of the energies at the scission point is 
shown in Fig. [SJ Due to the smaller Coulomb repulsion 
energy for more elongated shapes the scission point en- 
ergy calculated with the profile function © is by (2 -j- 5) 
MeV lower than that calculated with the profile function 

The expansion in multipole moments is an expansion 
in the complete set of orthogonal functions. At small 
deformation only few lower moments are important. At 
large deformation, especially at the scission point, one 
can not characterize the optimal shape of the surface 
by the quadrupole and octupole moments alone. Higher 
multipole moments should then also be taken into ac- 
count. 




0.4 0.6 
(M L -M R ) / (M L+ M R ) 

FIG. 5: Liquid drop energy ([5} at the maximal elongation 
R[ s 2 c * calculated with solutions of Eq.([8| (thick solid line), 
Eq. (|ll|l (thin solid line) and Eq. (|17[l (dash and dot lines) as 
function of the mass asymmetry 



The optimal shapes defined with R12, S constraints de- 
scribe well the separated or touching drops. For the 
shape with a neck one should try to define a constraint 
which would be an interpolation between Q2 , Q3 and 
Ri2,S constraints. Some hint how this can be achieved, 
one can get looking at the curvature of the surface calcu- 
lated with both constraints. At each point of the surface 
one can define the local curvature H(z), 



the profile function p{z) and the local curvature H{z), 
pp" = l + (pif-2pH(z)[l + (p'Y]i. (15) 

Comparing Eq. p3|) and Eq. (fTT]) one sees that the expres- 
sion in curly brackets in (J5J or pT|) is just twice the local 
curvature of the surface. 

The left-right asymmetric part of the curvature pro- 
portional to Signz and Ps{z/ ' zq) (the length of the drop 
along the z-axes is equal to 2zo) is shown in Fig. [5] The 
function P^^z/zq) which appears in Eq. (llO[) grows rapidly 
at the tips of the drop (at z ps zq). At the tips of the drop 
the surface of the heavy fragment becomes flat, the sur- 
face of the light fragment becomes very deformed (elon- 
gated). This is in contradiction with the expectation that 
due to Coulomb repulsion the distant parts of the drop 
should be close to spheroids. 

The spherical shape has a constant curvature. I.e. for 
the distant part of the drop the Sign z constraint is more 
meaningful than Q3. It is also clear that at small z (in 
the neck region) the curvature should change smoothly 
between the asymptotic values on the very left and on 
the very right. 
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FIG. 6: Illustration of the constraints S (solid), P^{z/zo) 
(dash) and smoothed Sign-function (thick solid line) 
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1 

R2 



(13) 



where R\ and R2 are the local principal radii of curva- 
ture. In the case of axially symmetric shapes the radii R\ 
and i?2 can be expressed in terms of the profile function 



i? 1 = P (z)Vi + (p') 2 , R2 = -[i + (p') 2 ]Vp", (14) 

Inserting (jT4"]) into (|13l) and solving this equation with 
respect to p"(z) one gets the following relation between 



These requirements can be fulfilled if instead of Sign z 
one would introduce a smoothed Sign - function, say by 
the replacement 



\z\ =^> Vz 2 + (Az) 2 , Signz =^ z/^/z 2 + (Az) 2 (16) 

and use these smoothed quantities as the constraints for 
the elongation and mass asymmetry. In (|16[) we have 
introduced also a smoother \z\ - function which appears 
in the definition ^ of R12 constraint. 
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IV. THE SCISSION SHAPES 



The replacement (|16p contains an additional parameter 
- the smoothing width Az. In principle, one can consider 
it as an additional collective parameter which has to be 
taken into account in the dynamical calculations. In the 
quasi-static limit one could expect that the value of Az 
is close the curvature radius i?2 in the neck region. 

Replacing in flTTT) the z and Signz by the smoothed 
quantities (|16p one gets the following equation for p(z) 



PP " = 1 + (p') 2 -p{X 1 V + A 2V /(z-z n )2 + (Az)2 + 



A, 



\/(z- Z n ) 2 + {Az)'- 



+ 10x LD ^ s (z) [I + (p') 2 }%17) 



In Figs. [JJ [5] we show the optimal shape at the scission 
point calculated with the constraint (I16I) for two very 
different values of Az, Az = 0.25i?o and Az = 2.0i?o- 
The first is approximately equal to the neck radius, the 
second - to the half-length of the drop in z-direction at 
the scission point. 



Q. 




1 
(z-z )/R n 

v cm' 

FIG. 7: The solutions of Eq. (JTTJ) (Az = 0.25R ) at the max- 
imal elongation ii^ f° r few values of the mass asymmetry 
5 = 0,0.1,. ..0.6 



The Figs. H] are similar to Figs. El SI The shapes 
calculated with Az = 0.25i?o are rather close the shapes 
calculated with Ri2,S constraints (in the limit Az — > 
the profile functions shown in Fig. [3] and Fig. [JJ coin- 
cide). The energies calculated with Az = 0.25i?o and 
i?i2, S constraints are also very close to each another, see 
Fig. With growing Az the scission shapes are get- 
ting shorter like those calculated with the QiiQj, con- 
straints. However the tips of the shapes calculated with 
Az = 2.0i?o are more "spherical" as compared with those 
calculated with Q2,Qs constraints. The energies of the 
shapes calculated with Az = 2.0i?o are by (1 -j- 3) MeV 
lower as compared with the ones calculated with Q2 , Q3 
constraints. 



Taking into account the results shown in Fig. [5j the 
shapes calculated with Az — 0.25i?o seem more prefer- 
able as the scission shapes. Besides, the total kinetic en- 
ergy of the fission fragments calculated with Az = 0.25i?o 
is in better agreement with the experimental data as 
compared with the one calculated for Az = 2.0i?o (see 
the Fig. Q3] below). So, in calculations below the shapes 
shown in Fig. [7] will be used as the scission shapes. 
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FIG. 8: The same as in Fig. ffj calculated for Az = 2.0i?o and 
6 = 0,0.1, ...0.8 

The description of the fission process requires a so- 
lution of the dynamical problem. The potential energy 
surface is important, but only one ingredient of the dy- 
namical description. The inertia, friction and diffusion 
tensors are also equally important. Still, having only the 
potential energy surface at ones disposal, one could try 
to estimate some observable of the fission process. 

Keeping in mind that the fission process is slow one 
could assume that during the fission process the state of 
the fissioning nucleus is close to thermal equilibrium, i.e. 
in the quasistatic limit the points {qi} on the deformation 
energy surface are populated with the probability given 
by the Boltzman factor, 



(18) 



Here T co u is the temperature, and Eq is a constant which 
is not important for what follows. 

Then, the normalized mass distribution of the fission 
fragments will be defined by the deformation energy at 
maximal deformation .R^ considered as a function of 
the mass asymmetry 5% for each pair of fragments, 



Y 



e -E def (R!& c, ->,S i )/T eoll 



(19) 



where E de f{R^ > 2 % '^ , ^i) is the liquid drop energy ([5]) plus 
the shell correction E s h e ii- 
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To calculate the shell correction energy E s h e ii we have 
approximated the shapes shown in Fig. [7] by the Cassini 
ovaloids with three deformation parameters a, ax, a^, see 
[H, and for the shape given in terms of Cassini ovaloids 
calculated the single-particle energies and the shell cor- 
rection by the code [3(. The liquid drop energy and the 
total energy including the shell correction for the nucleus 
236 U are shown in Fig. |H1 
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FIG. 9: The liquid drop (dash) and the total E to t = Eldm + 
E s heii deformation energy of 236 U along the scission line (at 
maximal elongation i?^"') calculated with the solutions of 
Eq.lQTJ for Az = 0.257? 



considered as a confirmation that the scission point shape 
and energy are calculated correctly. 

The distribution (fT5|) is a basic assumption of the 
scission-point model suggested by [14[ and developed fur- 
ther in |15l - [r7| . see also In this model the scision 
point configuration consists of two coaxial spheroids with 
tip-to-tip distance d and quadrupole deformation param- 
eters Pl and Ph- There are two temperatures in this 
model one for the population of the single-particle levels, 
Ti n t, and one for the collective degrees of freedom, T^oii- 
The parameters d, T int and T co u were fitted in [14j in 
order to reproduce the experimental data. The T co u was 
found to be close to 1 Mev. In the calculations shown in 
Fig. [TO] we used the same value T co u=\ MeV. 

Note, that within the optimal shape approach the 
shape configuration is not fitted to the experimental re- 
sults but is defined unambiguously from the minimal en- 
ergy condition. 



V. THE SEPARATED FRAGMENTS 

Within the optimal-shape method one can also find 
the optimal shape of separated fragments. For this aim 
one solves equation @ with the initial conditions that 
correspond to two spherical fragments at large enough 
distance R12 from each other. Making the distance R12 
smaller and smaller, one can find out how the shape of 
the fragments changes with the distance between their 
centers of mass. The results of numerical calculations 




7 / R 

fragment mass 



FIG. 10: The yield {I9j) of the fission fragments of 236 U cal- 
culated with the deformation energy shown in Fig. [5] 

The total energy has a minimum at the fragment mass 
equal to 95 or 141. Consequently the peaks of the mass 
distribution of the fission fragment are located at these 
fragments masses, see Fig.[SJ The position of the peaks of 
the mass distribution is in good agreement with the well 
known experimental results [13| . This agreement can be 



FIG. 11: The solutions of Eq. ((3]) corresponding to different 
values of the Lagrange multiplier A2 which fixes the deforma- 
tion R12 

for the symmetric splitting of 236 U are shown in Fig. [XT] 
The shape of the separated fragments is very close to 
oblate ellipsoids. The octupole deformation is very small. 
Its contribution to the deformation energy at the touch- 
ing point of two lls Pd nuclei is of the order 0.5 Mev 
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only. The energy of separated fragments is shown by the 
dashed line in Fig. [T^l The lower solid curve and the 
dashed line correspond respectively to the bottom of the 
fission and the fusion valleys and the upper solid curve 
to the ridge between the fusion and fission valleys. 

The kinetic energy of the fission fragments is the ki- 
netic energy gained by the fragments due to the Coulomb 
repulsion after separation plus the prescission kinetic en- 
ergy. Within the quasi-static picture one can calculate 
only the energy of the Coulomb repulsion "immediately 
after scission" . At present it is not so clear how the scis- 
sion process proceeds. For slow collective motion it is 
natural to assume that during the neck rupture the elon- 
gation (the distance between centers of mass of left and 
right parts of nucleus) does not change, like it is shown 
by arrow in Fig. 1121 The corresponding profile functions 

at i?i2 = R^'^ f° r the compact system and separated 
fragments are show in Fig. 1131 
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FIG. 12: Liquid-drop deformation energy ([5]) of 236 [/ as a 
function of the parameter R12 for the mono-nucleus (solid) 
and separated fragments (dash). 

The Coulomb interaction energy of the fragments im- 
mediately after scission shown in Fig. [13] is easy to calcu- 
late. The only additional parameter (besides the fissility 
parameter xld) which appears in such calculation is the 
parameter tq of the nuclear radius, Rq = tqA 1 ^ . In the 
present work we used the value ro = 1.225 fm. 

The comparison of the Coulomb interaction energy of 
the fragments immediately after scission with the exper- 
imental value of the total kinetic energy for the nucleus 
236 U is shown in Fig.rjJJ The solid and dashed lines are 
calculated with Az = 0.25i?o (solid) and Az — 2.0i?o 
(dashed). One sees that the more elongated scission 
(Az = 0.25i?o) shapes are in somewhat better agree- 
ments with the experimental data. 

The agreement is of qualitative character only. For a 
more accurate description on should take into account 
the multimodal character of the fission of 236 U. The 
optimal shapes described above correspond to only one 
(standard) fission mode. 
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FIG. 13: Profile function, solution of Eq. (|17|) at the scission 
point R[ s 2 c1 for a mass-symmetric deformation of 236 U (solid 
line) , and profile function ([3| of two separate fragments with a 
distance between their centers of mass equal to ii^ (dashed 
lines) 
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FIG. 14: Coulomb repulsion energy of two separated frag- 
ments at a distance R12 = R^*'\ where R[ a 2 cl is the max- 
imal elongation calculated with Eq. (|17[) with Az = 0.25i?o 
(solid line) and Az — 2.0i?o (dashed line). The experimen- 
tal values of the total kinetic energy are taken from [T^-[2ll] 
The dash-dot and dash-dot-dot lines show Coulomb repulsion 
energy calculated for super-long and super-short shapes (see 
Fig.Q]| 



VI. SUPER-LONG AND SUPER-SHORT 
SCISSION SHAPES 

The optimal shapes discussed in Sections 3-5 have the 
two degrees of freedom - elongation and the mass asym- 
metry. The neck radius for the given elongation and the 
mass asymmetry attains the " most favored" value which 
results from the minimum of the potential energy condi- 
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tion. In the dynamical calculations of the hssion process 
the neck radius is often considered as an independent 
collective variable which can deviate from the one corre- 
sponding to the bottom of the potential energy surface. 
Thus, it makes sense to incorporate in the optimal shapes 
procedure the neck radius as another independent degree 
of freedom. 

In order to include one additional degree of freedom in 
the optimal shapes procedure one should add another 
constraint fixing the neck radius. Usually, in various 
shape parameterizations, the neck radius is regulated by 
the parameter of the hexadecapole deformation. Using 
X4Q4 as an additional constraint, allows, indeed, to vary 
somewhat the neck radius of the drop. However, at large 
value of A4 the A4Q4 constraint results in very peculiar 
shapes. 

Another possibility to vary the neck radius is to fix the 
amount of matter in the neck region by introducing the 
constraining function f± of the type 



terpretation of the experimental data, see for example 
23]. 



1 

V 



dVp 2 (z) exp 



(20) 



For simplicity we assume here that Az has the same 
meaning and value as used in Sections 3-4. 

The effect of A4Q4 on the optimal shapes is demon- 
strated in Fig. [T5] Indeed, varying A4 (keeping A2 fixed) 
allows to change the neck of the drop in a rather broad 
region. 
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FIG. 15: Effect of a A4/4 restriction (|20|l on the optimal 
shapes: the heavy solid line shows the profile function cal- 
culated with A4 = 

The introduction of the neck degree of freedom has an 
important consequence for the scission shape. Depending 
on the neck radius, the scission shapes become more elon- 
gated or shorter, see Fig. [T6l Thus, it turns out possible 
to introduce the so called [22| super-long or super-short 
scission shapes which represent the possibility of the ex- 
istence of few fission modes and are exploited by the in- 



Q. 




z/R 

FIG. 16: The profile functions p(z) at the scission point calcu- 
lated with A4 = (solid), A4 = 0.5 (dash-dot) and A4 = —0.5 
(dash-dot-dot). 

The Coulomb interaction energy "immediately after 
scission" for the super-long or super-short shapes shown 
in Fig. [T5] is plotted in Fig. [TJ] by dash-dot and dash- 
dot-dot lines. Qualitatively these results are very close 
to the contribution from three fission modes 24] shown 
in Fig. 13 of [IE)]. For a more precise estimate of the 
contribution from the super-long or super-short scission 
shapes, full dynamical calculations (with the account of 
shell effects) are required. 



VII. SUMMARY AND CONCLUSIONS 

The optimal-shape approach is put into practice by 
the construction of the constraint on the mass asymme- 
try which is an interpolation between the constraint on 
quadrupole and octupole moments (which is quite suc- 
cessful at small deformations) and the constraint on the 
distance between the centers of mass of the future fis- 
sion fragments and the difference of their masses (which 
is well defined for the shape with a neck or separated 
fragments). The use of this new constraint allows to de- 
fine the scission point shapes in broad region of the mass 
asymmetries. 

It is shown that the optimal-shape procedure can be 
further extended by incorporating the neck degree of free- 
dom. The introduction of the neck degree of freedom 
leads to the fission valleys, the existence of which follows 
from the analysis of experimental data. 

The account of shell effects on the potential energy 
surface of the optimal drops will be the subject of future 
studies. 
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